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ABSTRACT 

Two general types of numerical models for predicting microwave emission from soils are com- 
pared -coherent and noncoherent. In the former, radiation in the soil is treated coherently, and 
the boundary conditions on the electric Fields across the layer boundaries are used to calculate the 
radiation intensity. In the latter, the radiation is assumed to be noncoherent, and the intensities 
of the radiation are considered directly. The results from the two approaches may be different be- 
cause of the effects of interference, which can cause the transmitted intensity at the surface 
(i.e., emissivity) to be sometimes higher and sometimes lower for the coherent case than for the 
noncoherent case, depending on the relative phases of the reflected fields from the lower layers. 
This coupling between soil layers in the coherent models leads to greater soil moisture sampling 
depths observed with this type of model, and is the major difference that is found between the two 
types of models. In noncoherent models, the emissivity is determined by the dielectric contract at 
the air/soil interface. The subsequent differences in the results are functions of both the frequency 


of the radiation being considered and the steepness of the moisture gradient near the surface. The 


calculations were performed at frequencies of 1.4 and 19.4 GHz and for two sets of soil profiles, 
little difference was observed between the models at 19.4 GHz; and only at the lower frequency 
were differences apparent because of the greater soil moisture sampling depth at this frequency. 



A COMPARISON OF RADIATIVE TRANSFER MODELS FOR PREDICTING 


MICROWAVE EMISSION FROM SOILS 


T. J. SCHMUGGE 
Hydrological Sciences Branch 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 

and 

B. J. CHOUDHURY 
Computer Sciences Corporation 
Silver Spring, Maryland 


INTRODUCTION 

Soil moisture information is important in a variety of disciplines (e.g., hydrology, agriculture, 
and meteorology). Attempts to monitor world food supply, predict watershed runoff, and model 
boundary layer heat and moisture conditions are but a few of many specific areas in which surface 
and subsurface soil moisture information is needed. To this end, both active and passive micro- 
wave remote sensing techniques are being studied (Ulaby et al., 1978 and 1979; Schmugge, 1978) 
to provide efficient and cost-effective means of estimating average moisture content in wide areas. 

A key parameter in understanding and interpreting remotely sensed data is the dielectric constant 
of the target medium, which has a large effect on the reflective and emissive properties of the 
surface. Microwave frequencies are ideally suited for soil moisture remote sensing because at these 
frequencies the dielectric constant changes rapidly with moisture content. In fact, radiance meas- 
urements obtained by ground-based, airborne, and satellite-borne radiometers have indicated 
emissivities from as low as 0.6 when wet to 0.95 when dry. Microwave radiative transfer models 


have been developed to provide the theoretical basis for this relationship. 


Three recent papers have presented theoretical models for microwave emission from soils (tyjoku 
and Kong, 1977; Wilheit, 1 978; and Burke, et al., 1979). These models considered die emission 
from the soil lor a range of moisture and temperature profiles and studied die effects of varia- 
tions in these subsurface properties on the emission from die surface. The purpose of this paper is 
to compare these models in terms of their fundamental approaches and their remits for a set of 
standard soil profiles. The models will be compared at the microwave frequencies of 1.4 and 19.4 
GHz (wavelengths of 2 1 and l.SS cm). 


For homogeneous media (i.e„ those with constant moisture and temperature profiles), die emissivity 
(e) for a smooth surface can be calculated from the soil’s dielectric properties using die Fresnel 
equations for surface reflectivity (r). For perpendicular incidence, these equations reduce to: 

IvT-il’ 


e * 1 - r ■ 1 - 


v7+ 1 


( 1 ) 


where e is the dielectric constant of the soil. 


Difficulties arise when both the temperature and moisture content of the soil vary with depth in 
the soil. The basic approach for solving this problem is to break the sod volume into thin layers and 
then numerically sum the contributions of each layer. The various ways of treating this numerical 
summation through the soil layers form the basis of major differences between the particular 
radiative transfer models. 

There are two alternate approaches to studying microwave emission from soils based on the assumed 
characteristics of radiation within the soil (i.e., whether it is treated coherently or noncoherently). 
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If radiation within the soil is assumed to be coherent, its intensity must be obtained by calculating 
the electromagnetic field vectors from a solution to Maxwell’s equations. However, if radiation 
within the soil is assumed to be incoherent, the intensity of radiation can be considered directly. 
The type of assumption made will affect the performance of the models since intensities calculated 
for coherent and incoherent radiations may be different because of wave interference effects 
associated with coherent radiation. These interference effects cause the coherent intensity to be 
sometimes higher and other times lower than the incoherent intensity. Oscillations of the intensity 
as a function of frequency is an indication of coherency effects in the radiation. 

« v. 

The need for using a coherent model to study thermal emission that is intrinsically a noncoherent 
process may be questioned. However, if a spatial filter (e.g., a dielectric slab) is introduced, the 
emerging radiation will be partially coherent, which should be manifested in interference phenomena. 


Indeed, radiometric observations of interference effects have been obtained for sandbox experi- 
ments in which a metal plate is buried by increasing thicknesses of sand (Blinn, et al., 1972) 
and for oil films on water (Hollinger and Mannella, 1973). In an analysis of this problem, Carver 
(1977) concluded that if the layer thicknesses were less than the coherence length of the emitted 
photons (L) given by 


L = c/(Af /l) 


( 2 ) 


where c is the velocity of light, Af is the bandwidth of the radiometer, and e is the dielectric con- 
stant of the medium, then coherence effects are important. At 1.42 GHz, Af is generally 27 MHz, 
the bandwidth of the radio astronomy band at this frequency, and e of soil will range from 4 to 25, 
yielding a value of L between 5 and 2 meters. This range is much greater than that of any of the 
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layer thicknesses used In the models considered here and thus coherent effects may be important 
The thicknesses were chosen to be much less than the distance over which there was a substantial 
change in soil properties. These major changes in soil moisture content occur in distances much 
less than I.. For these reasons, it is necessary to consider coherent models. As will be noted, 
coupling between reflectivity at the surface and tiie subsurface dielectric properties occurs hi the 
coherent models but not in the noncoherent models and leads to the major differences between 
the two approaches. 

All three models assume that scattering by pebble* or other discontinuities in the soil is insig- 
nificant. In a study of the microwave emission from a scattering medium, England (1975) found 
that scattering effects will be significant when the scattering albedo (co 0 ) is greater than 0. 1 , where 

co c , the ratio of the scattering loss to the total loss, is given by: 

No 

w » 

° No + 2 ft (3) 

where a is the scattering cross section of the pebble,)) is the loss in the medium, and N is the 
number of pebbles per unit volume. Assuming appropriate values for the dielectric properties of 
pebbles and wet soils and assuming a fractional pebble volume of 1 percent, it can be found from 
England’s equations that volume scattering will be sipuficait only if d/X > 1/7, where d is the 
pebble size and X is the wavelength (0,5-cm pebbles in the 3-cm or X-Band wavelength region). 

Even if the size or number of particles becomes larger, volume scattering would have little effect 
if X > 10 cm and will therefore be neglected in this paper. 

In the next section, the three models to be considered will be briefly described. The results obtained 
using these models on the soil profiles published in the Njoku and Kong paper are presented and 
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compared in Section 3. Further calculations were performed using the Burke, Wilheit, and a 
simpler radiative transfer model for some actual profiles But possessed steep moisture gradients 
near die surface. These results are also presented in Section 3. 


COMPARISON OF THE MODELS 
Noncoherent Models 

The microwave emission model developed by Burke, et al. (1979), hereafter referred to as the Burke 
model, is based on die assumption of incoherent radiation. Radiation intensity is calculated from 
solution of the radiative transfer equation: 

dT r 

— - * 7<2)T p + y(*)T(2) (4) 

where y is the absorption in the layer, T p is radiation intensity expressed in terms of a temperature, 
T is the physical temperature of die soil, z is the path length, and subscript p indicates die polariza- 
tion of the radiation. 


Soil is treated as a layered dielectric; each layer is homogeneous and of arbitrary thickness. Thick- 
ness (Azj), absorption coefficient ( 7 j), and temperature (Tj) of each soil layer are specified. 

The radiative transfer equation for each layer is then solved. The radiative transfer equation in the 
tint layer may be written as: 


dT p 

— — - . T + T 

d(7,z) 9 1 


(5) 


This equation can be integrated from a point just below die surface to a point just above the inter- 
face between the first and second layers. Because the dielectric properties are assumed to be con- 
stant across the layer. 


T p (l ) - T, (1 - e' 1, i A *i) ♦ T p (2*)e-»i A, i 


( 6 ) 
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The argument (N±) implies that the measurement is made above (plus) or below (minus) the Nth 
interface. The first term on the right-hand side accounts for radiation emitted within die first layer 
which comes directly to the surface. Hie second tom describes upwelling radiation at the bottom 
of the first layer. This in turn has two components: ( 1 ) radiation emitted in die first layer and re- 
flected at the interface between the first and second layers; mid (2) radiation transmitted from 
lower layers. 

T p (2 + ) - T, (1 - e^ A *i) ♦ T p (2'Xl • R*) (7) 

R p2 is the absolute value squared of the Fresnel coefficient for the p polarization. The radiation 
field just above the surface is the value just below multiplied by the transmittance ( 1 - Rp t ) 

yi + ) Mi-yvn (8) 

- 0 - R pI )|t,( 1 - e' T » 4 *t)(l + Rpje-V*,) ♦ (1 - R p3 )Tp(2)e>i A *i( 

The radiative transfer equation can be integrated again to calculate T p (2*). Repeating the procedure 
for N layers gives the brightness temperature as: 



(9) 


Beginning with the deepest layer of the soil, the intensity emerging from each layer is calculated to 
obtain the observed microwave intensity. Results from this model were compare* with actual ob- 
servations at the 2.8- and 21 -cm wavelengths (Burke, et al., 1979), with layer thickness correspond- 
ing to observed soil moisture values in sampling depths of 0 to 1, 1 to 2, 2 to 5, 5 to 9, and 9 to 
IS cm from the surface. Some of the conclusions from the study were that: (1) the emitted 
intensity at both wavelengths correlates best with the near-surface moisture;(2) the slope of the 
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intensity/moisture curves decrease in going from day to dawn; and (3) i nc reased neanurface 
moisture at dawn is characterized by increased polarization. 


Coherent Models 

Microwave emission models developed by WUheit (1978) and by Njoku and Kong (1977), hereafter 


referred to as the Njoku model, are based on the assumption of coherent radiation. In Wilheit's 
model, soil is treated as a layered dielectric. Solutions of Maxwell's equations and the boundary 
condition at the interfaces are used to calculate the electric field in each layer. These electric-field 
values are used to calculate the energy fluxes and thus the fractional absorption in each layer 


f f 




( 10 ) 


where S£j is the electromagnetic field energy (i.e., Poynting vector) entering the jth layer at the 
0-1 )th interface, Sj* is the energy for the 0+1 )th layer at the jth interface, and SJ is the energy 
incident on the first interface. The superscript (p) designates the polarization of radiation. These 
energy values are calculated using Poynting’s theorem for electric fields. If Tj is the temper- 
ature of the jth layer, under thermodynamic equilibrium, the layer radiates energy equal to the 
product of the fractional absorption (fp and the temperature (Tj). The brightness temperature 
is given by 


n ■ £ t 


’ (id 

The conservation of energy at the air/soil interface determines the reflectivity of the soil R p as: 


R p * l 



( 12 ) 


7 


Computations with the model have indicated that radiation from the soil is characterized by 
two sampling depths: reflective and thermal (Wilheit, 1978). The reflectivity is characterized by 
changes in the real part of the index of refraction over a sampling depth: $ r 0.1 X, where X is the 
wavelength in the medium. The thermal sampling depth is determined by the imaginary part of the 
index in the medium and is given by: 



(13) 


where Xj is the depth of the ith layer. For a uniform dielectric, this reduces to 

» X 

T * 4jt lm(n) 

For a low-loss dry soil, 6 T will be an order of magnitude larger than 6 t , whereas for a wet soil, it 
will be only slightly larger. 


Another parameter of interest is the average soil temperature over this thermal sampling depth, 
which is referred to here as the effective radiating temperature of the soil and is given by: 


T|f, 


T eff “ 


(15) 


where T is the physical temperature of the ith layer. The ratio of T B to T eff is an effective emis- 
sivity for the soil. 


Recently, this model, incorporating the effect of surface roughness on the transmission coefficient, 
was used to explain observed emission intensities for 1 .55- and 21 -cm radiations (Choudhury, et a!., 
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1979). Although the agreement was generally satisfactory, in some cases the wave interference 
effect may have ' ■■id to the prediction of higher intensities than those observed. 

The Njoku model is based on the coherent radiative transfer formulation of Stogryn (1970) for a 
thermally and dielectrically nonhomofeneous medium. In tills formulation, the fraction of energy 
absorbed at different depths within the medium is obtained from the solution of a differential 
equation with a flux conservation boundary condition at the air/soil interface. This method of 
calculation is the same as the one used for calculating the reflectivity of a nonhomofeneous dielec- 
tric, By integrating the product of fractional absorption and black body emission intensity, the 
observed intensities are obtained. For a set of simulated moisture and temperature profiles, Njoku 
and Kong' (1977) have given illustrative results for wavelength dependence of the emitted intensity. 

A continuous nonhomofeneous dielectric can be well approximated by a layered dielectric, using 
appropriate choice of layer thicknesses. It may thus be expected that Wilheit and Njoku's models 
will closely agree with one another because both are based on the assumption of coherent radiation. 
However, the method of calculating the observed intensities are quite different in the two models. 
Physically and numerically, Wilheit's model is simpler and more efficient The following section 
presents a numerical comparison of these models. 

In the algorithm used for Wilheit's model, the soil was divided into 100 layers whose thicknesses 
varied with soil depth. Beginning at the surface, there are 1 0 layers having each of the thicknesses 
0. 1 , 0.2, 0J, and 0.4 cm, followed by 20 layers O.S-cm thick and 40 layers 1 ,0-cm thick for a total 
soil depth of 60 cm. With this soil grid, moisture and temperature values are obtained from 
the analytical expressions given by Njoku and Kong, and the dielectric constants are interpolated 
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from the values presented by them. The brightness temperature and effective radiating temperature 
are calculated using the Fortran routines provided in Wilheit’s paper. 

In the Burke model, there are also 100 soil layers whore thicknesses varied exponentially with 
depth. The surface layer thickness is about 0.003 cm, increasing to 1 cm at a depth of 9 cm and 
5 cm at a depth of 40 cm. The total soil depth in this model is 100 cm. 


QUANTITATIVE COMPARISONS 

To compare the three models quantitatively, values of T B were calculated for both the Wilheit and 
the Burke models, using the soil profiles given in the Njoku ard Kong paper. The values of T B 
obtained in this way were compared to those interpolated from the figures in the Njoku and Kong 
paper. The moisture (p(z))and temperature (T(z)) profiles used are given by the following equations: 

p(z) ■ P + Ap — - — for 0 < z < d 

‘ [e* d - 1] (16j 

» p(d) for z >u 


T(z) - T 0 + AT, e“»i‘ ♦ AT 2 «' Y 2* 


(17) 


Table 1 lists the values of the parameters used in these equations, and Figure 1 graphically presents 
the resulting moisture and temperature profiles. These profiles cover a range that is representative 
of those that could realistically be encountered in nature. 


The soil dielectric constants used were those given in the Njoku and Kong paper. Therefore, the 
results obtained using the Burke and Wilheit models should be comparable to those given by Njoku 
and Kong. The calculations were performed at 1.4 and 19.35 GHz (21 and 1.55 cm) frequencies. 
Comparisons at both frequencies appear in Tables 2a through f. 
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Table 1 

Parameters of the Moisture and Temperature Profile* 
(from Njoku and Kong, 1977) 


Moisture Profile 


Ap 


d,m 


30.0 -5.0 10.0 0.5 

15.0 10.0 50.0 0.5 

18.0 20.0 0.5 

18.0 5.0 0.5 

10.0 30.0 0.5 

2.0 20.0 -10.0. 0.5 


Temperature Profile 


300.0 


AT, 

AT a 

15.0 

0.0 

-15.0 

1 

-15.0 

20.0 

0.0 

0.0 

-30.0 

0.0 

15.0 

0.0 


At both frequencies, the agreement among the three models is good; die maximum difference is less 
than 4 K at 1.4 GHz and about 2 K at 19.4 GHz. There are, however, significant differences between 
the results at the two frequencies. For the wetter profiles (numbers 1 and 2), the 1.4 GHz has lower 
T b ’s because of the larger dielectric constant for wet soils at this frequency. Also die lower fre- 
quency, longer wavelength, has a greater temperature sampling depth. This tendency can be observed 
by comparing the differences between temperature profi! : 1 and 2 for the two frequencies, which 
have a 3 K difference at the surface but are equal at about 35 cm. This is particularly noticeable for 
moisture profile 4, the driest, where there is a 5 K difference in response at 1 .4 GHz compared with 
22 K at 19.4 GHz. 


The weighting functions for each model were also determined; the results agreed with those presented 
in the Njoku and Kong paper and will not be presented here. Since the weighting functions describe 
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Figure 1. Moisture and temperature profiles representative of 

various soil-moisture conditions. Profiles correspond to 
parameters of Table 1 (from Njoku and Kong, 1977). 



Table 2 

Calculator Brightness Temperatures for the Njoku and Kong Soil Profdes 


(a) Temperature Profile #1 


Moisture 

Profile 


(b) Temperature Profile #2 


Moisture 

Profile 


1 .4-GHz Frequency 

Burke 

Wilheit 

Njoku 

191.9 

192.0 

193.0 

233.0 

230.5 

230.0 

267.3 

263.5 

264.0 

278.9 

280.5 

277.5 


19.4-GHz Frequency 
Burke I Wilheit I Njoku 


203.6 

242.3 

279.3 

289.7 


203.6 

241.9 

279.5 

289.6 


206.0 

242.0 

280.0 

289.0 


Moisture 

Profile 


Moisture 

Profile 


1 ,4-GHz Frequency 
Burke 1 Wilheit I Njoku 


184.5 

225.4 

260.4 

274.2 


184.6 

223.0 

256.7 
275.9 


185.0 

222.5 

257.0 

272.5 


19.4-GHz Frequency 

Burke 

Wilheit 

Njoku 

184.5 

184.6 

187.0 

219.9 

219.5 

219.0 

254.5 

254.7 

256.0 

267.1 

266.9 

267.0 


(c) Temperature Profile #3 


Moisture 

Profile 


1 ,4-GHz Frequency 

Burke 

Wilheit 

Njoku 

187.7 

187.8 

188.0 

228.4 

226.0 

226.0 

262.8 

259.1 

260.0 

275.7 

277 3 

274.0 


(d) Temperature Profile #4 

Moisture 1.4-GHz Frequency 

Profile Burke | Wilheit I Njoku 


188.2 

229.2 

263.8 

276.5 


188.3 

226.8 

260.1 

278.2 


189.0 

227.0 

261.0 
275.0 


Moisture 

Profile 


!'• 1-GHz Frequency 

Burke 

Wilheit 

Njoku 

197.1 

197.1 

199.0 

234.5 

234.1 

234.0 

270.0 

270.2 

271.0 

279.8 

279.7 

278.0 


Moisture 19.4-GHz Frequency 
Profile Burke I Wilheit I Njoku 


194.1 

231.1 

266.9 
278.4 


194.1 

230.7 

267.1 

278.3 


196.0 

230.0 

268.0 

276.0 





Table 2 (continued) 


(e) Temperature Profile #5 


Moisture 

Profile 


Burke 

Wilheit 

Njoku 

1 

185.9 

186.0 

186.0 

2 

227.6 

225.2 

225.0 

3 

263.6 

260.0 

261.0 

4 

279.5 

281.2 

278.0 


Moisture 

Profile 

19.4-GHz Frequency 

Burke 

Wilheit 

Njoku 

1 

184.5 

184.5 

187.0 

2 

219.8 

219.5 

219.0 

3 

254.3 

254.5 

256.0 

4 

266.6 

266.5 

267.0 


(f) Temperature Profile #6 


Moisture 

Profile 

1.443Hz Freauency 1 

Burke 

Wilheit 

Njoku 

5 

6 

268.8 

266.4 

266.2 

270.4 

266.0 

268.5 


Moisture 

19.4-GHz Frequency 

Profile 

Burke 


Njoku 

5 

281.5 

281.7 

280.0 

6 

279.6 

279.6 

278.0 


the depth distribution of the sources of the radiation, this agreement indicates that all the three 
models yield the same effective radiating temperature and that any differences should arise from 
their effective emissivities. This can best be seen using the results for constant temperature profile 
number 4 in Table 2d. The emissivity can be obtained directly by dividing the calculated T B by 
the soil temperature, 300 K. These results are given in Table 3. 


The biggest difference between the coherent and noncoherent models appears for moisture 
profiles with the steepest gradient near the surface (i.e., for profiles 2 and 3 in Table 3 emissivity 
results from the noncoherent model are greater by 0.01 at 1.4 GHz, with no difference between 
the models at 19.4 GHz) because the emissivity calculated in the Burke model corresponds to the 
Fresnel transmission coefficient calculated for the surface moisture value. Since it is predominantly 
influenced by surface moisture, the Burke model emissivity will increase with decreasing surface 
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Table 3 

Emissivities Determined from the Brightness Temperature 
Results for Profile 4 (T(z) ■ 300 K) 


Moisture 

1.4-GHz 

19.4-GHz 

Profile 

Burke 

Wilheit 

Njoku 

Burke 

Wilheit 

Njoku 

1 

0.627 

0.628 

0.630 

0.647 

0.647 

0.653 

2 

0.764 

0.756 

0.757 

0.770 

0.769 

0.767 

3 

0.879 

0.867 

0.870 

0.890 

0.890 

0.893 

4 

0.922 

0.927 

0.917 

0.928 

0.928 

0.920 


moisture, irrespective of the nature of the actual soil moisture profile. In coherent models, however, 
the emissivity is dependent on the moisture profile just below the surface because the phase inform- 
ation introduces a connection between the surface and the soil layers beneath it. Thus, the 
coherent models will depend on the average moisture content of the soil moisture sampling depth 
(which Wilheit estimates to be about one-tenth of the wavelength in the medium). Therefore, the 
agreement between the coherent and noncoherent models would be better if the thickness of the 
surface layer could be chosen to be equal to this sampling depth. 

To study the effect of these steep moisture gradients in more detail, calculations were performed 
for a series of soil moisture profiles measured at the U. S. Water Conservation Laboratory in 
Phoenix (Jackson, 1973) as a field dried out after an irrigation in March 1971. The profiles 
presented in Figure 2 are midday profiles that had the steepest moisture gradients near the surface. 
In addition to the Burke and Wilheit models, a variant of the Burke model was also investigated. In 
this model, the Fresnel reflectivity at all layer interfaces except for the soil/air interface was set to 
zero (i.e., it was assumed that all the soil interfaces were transparent). The rationale behind this 
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SOU DEPTH (Cm) 
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SOIL DEPTH (Cm) 
(b) MOISTURE 


Figure 2. Soil-moisture and temperature profiles observed at the U. S. Water 
Conservation Laboratory. The curves are labeled by the number of 
days after irrigation. 
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DIELECTRIC CONSTANT 



SOIL MOISTURE (WEIGHT-%) 


Figure 3. Real and imaginary pails of the dielectric constant for soils with textures 1 

similar to those of the soil at the U. S. Water Conservation Laboratory. The j 

data are from laboratory measurements at the frequencies of 1 .4 (X ■ 21 cm) 1 

and 19.4 (X= 1.55 cm) GHz. j 
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model is that, if the layer thicknesses are chosen so that the dielectric constants of two adjacent 
layers are nearly equal, the interface reflectivity will be negligible. For this model to be valid, the 

t 

layer thicknesses should be reduced as the moisture profile gets steeper. The brightness temperatuie 
in this model is given by: 


TP 


(1 - R p ) J T(a) 7 (*)«p ^- £ y(*')dz^ dz 


(18) 


(19) 


Ml -tyT* 

Note that, as in the Burke model, the first factor on die right side is the Fresnel transmissivity of 
die polarization (p) calculated for the surface moisture value, and die second factor is die direct 
integration of the radiative transfer equation (1) for the intensity just below the soil surface which 
yields the effective radiating temperature (T eff ) of the medium. The layer thicknesses used for the 
Burke model were also used for this model. This model was presented in equation 1 0 of the Njoku 
and Kong paper and will be called the radiative transfer model. The results of Njoku and Kong 
indicate that there is good agreement between the radiative transfer model and their coherent model 
for frequencies above about S GHz but that differences arise for lower frequencies. This difference 
will be illustrated further by calculations made using the seal profiles given in Figure 2. The dielec- 
tric constants used with these profiles appear in Figure 3. These data are for soils with textures 
similar to that of the soil at the U. S. Water Conservation Laboratory. 




Table 4 gives the results for the three models (Burke, Wilheit, and radiative transfer) in terms of 
both T B and emissivity (e = T B /T eff ). Note the excellent agreement between the Burke and 
radiative transfer models for all profiles. There are only small differences between the Burke and 
Wilheit results for the 19.4-GHz calculations, but there are substantial differences in tire 1.4-GHz 
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Table 4 

Brightness Temperatures Calculated for the Profiles in Figure 2 at 1.4 and 19,4 GHz.* 


Profile 

1.4 GHz 1 

! 19.4 GHz 

Radiative Transfer 

Burke 

Wilheit 

Radiative Transfer 

Burke 

Wilheit 

3 

177.3 

177.5 

173.2 

191.0 

191.2 

195.6 

4 

203.9 

204.1 

202.7 

212.6 

212.9 

217.1 

5 

250.7 

249.2 

229.4 

272.0 

271.6 

275.5 

6 

262.0 

260.4 

242.1 

279.9 

279.2 

281.2 

8 

270.3 

269.0 

258.4 

285.0 

284.4 

286.8 

12 

271.6 

270.5 

268.3 

284.7 

284.3 

287.0 

16 

273.6 

272.6 

273.0 

287.5 

287.2 

288.7 

23 

273.9 

272.9 

276.1 

290.4 

290.4 

291.1 

37 

278.9 

278.0 

279.1 

295.7 

295.7 

296.3 


•The profiles are identified by the number of days after irrigation. 

calculations-up to 20 K for the profiles with the steepest gradients (i.e., profiles 5, 6, and 8). The 
good agreement of the 19.4-GHz results indicate that the layer thicknesses used for the noncoherent 
models are approximately the same as the sampling depth. The differences in the 1.4-GHz results 
indicate that the sampling depth is greater than the first layer thickness used in the noncoherent 
models. 


The 1 K differences betv 'een the radiative transfer and Burke models indicates that, although the 
reflectivity is small between the soil layers, there is some minor effect. However, as shown here at 
1.4 GHz, the effect would generally be negligible. 

The emissivities for these profiles can be determined by using the effective radiating temperature 
as defined by equations IS and 19. Because both the Burke and Wilheit models were found to have 
the same weighting functions, the effective temperature could be calculated by using either model, 
and in this case the results from the Wilheit model were used. Table 5 presents the values of T eff 
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plus the subsequently derived emissivities. The measured surface (depth ■ 1 mm) temperatures 
are listed for comparison. 

The effect of the greater temperature sampling depths at 1.4 GHz is apparent in the values of T aff . 
At 1.4 GHz, there is a 10 K increase from wet to dry (profiles 3 and 37) and at 19.4 GHz there is a 
21 K increase for the same 29 K increase in surface temperature. For the wet profile, the values of 
T tff at 19.4 GHz are close to the surface temperature, but, as the soil dries, the two temperatures 
diverge because the microwave thermal radiation originates at deeper depths in the soil where it is 
cooler for these profiles. Table S also gives the thermal sampling depths defined by equation 13 for 
these two frequencies. 

On the days with the steepest moisture gradients (i.e„ 5, 6, and 8), the largest differences between 
the models occur because, as noted previously, the Burke model uses the moisture content of the 
surface layer only to determine the emissivity. These differences cannot be accounted for by the 
fact that different layer thicknesses were used for our calculations with the two models. The 
moisture content in the 1-mm surface layer used for the Wilheit model is essentially the same as 
that for the 0.03-min layer used in the Burke model. In fact, it was the thinness of the layers used 
in the noncoherent models that afforded the excellent agreement between the radiative transfer and 
Burke models. The difference between the Burke and Wilheit models could be reduced if a thicker 
surface layer were used in the calculations with the noncoherent models. 

An estimate of this desired thickness can be obtained by comparing the emissivities obtained with 
the Wilheit model with those calculated using the equation for a uniform dielectric. Emissivities 
calculated using the Fresnel equations are given in Figure 4 based on the dielectric values presented 
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in Figure 3. For example, on day S the emissivity calculated with the Wilheit model is 0.789, which 
intersects the curve in Figure 4 at a moisture of 12.5 percent for a uniform soil. Table 6 lists the 
remits of this comparison for days 4, 5, 6, 8, and 12 and the average soil moistures in the 0* to 1-cm 
and the 0- to 2.5-cm layers of the soil. 

In Table 6, the Fresnel moisture is obtained from the intersection of the model calculated emissivity 
with the curve in Figure 4. For days 4, 5, and 6, there is good agreement of the Wilheit value with 
the average measured moisture in a 0- to 2.5-cm layer. As would be expected, the Burke value 
is drier than the 0- to 1-cm layer for all profiles. An interesting feature is that, for the drier 
days (8 and 12), the Wilheit value indicates a moisture as dry or drier than the 0- to 1-cm lav r. 

This could result from a positive interference effect increasing the effective transmission across the 
surface. Interference effects undoubtedly cause the higher emissivities calculated with the Wilheit 
model compared with Burke model for days 16, 23, and 37. This effect was also noted in some 
other calculations with the coherent model (Choudhury, et al., 1979). 


Table 6 

Comparison of Model and Fresnel Emissivities 


Day 

Wilheit 

Emissivity 

Fresnel 

Moisture 

Burke 

Emissivity 

Fresnel 

Moisture 

Soil Moisture (weight-%) 

0 - 1 cm 

0-2.5 cm 

n 

0.699 

15.7 

0.704 

15.3 

16.1 

16.5 

D 

0.789 

12.5 

0.857 

I 

11.0 

13.8 

H 

0.829 

10.8 

0.892 


8.4 

12.0 

8 

0.877 

6.0 

0.913 

1 

5.6 

9.9 

12 

0.910 

3.0 

0.918 

2.1 

4.4 

8.3 
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SOU M04STUM (WIMMT-%) 

Figure 4. Calculated values of the emissivity at 1.4 GHz using the dielectric constants 

from Figure 3 in Equation 1. 
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CONCLUSIONS 


Several models for estimating microwave emtaton from toils were compared-the htfoku, VO* 
heit, Burke, and radiative transfer models. On the basis of their fundamental approach, they cm be 
classified into noncoherent (Burke and radiative transfer) and coherent (Njoku and Wilheit) models; 
die former is ’' seed on die intensity and die latter mi die amplitude of the radiation (fold. Assoc- 
iated with this difference in fundamental approach, results obtained from these two types of models 
differed, especially in regard to their apparent soil moisture sampling depth. In the noncoherent 
models, emissivity is determined by the dielectric contrast at the air/soil interface; therefore, for a 
soil with a steep moisture gradient at the surface, this emissivity will depend on layer thicknesses 
used. Because the coherent models keep track of the phase of the fields in the soil, the value of the 
emissivity at the surface is coupled to the dielectric properties of the layers below die surface. 

Due to the dependence of the emissivity in die noncoherent models on the air-soil dielectric con- 
trast, these models do not predict accurate soil moisture sampling depths. The magnitude of 
the difference of this sampling depth between the coherent and noncoherent models wUI depend on 
the microwave frequency and the steepness of die soil -moisture gradient in the soil. Njoku and. 
Kong showed that, for their profiles, there was no difference for frequencies above 4 GHz. 

As noted previously, interference effects cm occur in the coherent models; however, these effects 
are not likely to be observed in nature because of the roughness of the surface and horizontal 
inhomogeneitics in the soU. Therefore, care mutt be exercised when interpreting the result! from 
the coherent models. 
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All of the models appeared to have the tame weighting functions for distributing the sources of the 
radiation, and, as a result, they all have the same thermal sampling depth and predicted die same 
effective radiating temperature. 
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